LABORATORI NAZIONALI Dl FRASCATI 

SIS-Pubblicazioni 



LNF-04/14 (P) 

28 luglio 2004 

hep-ex 



NEW CHARM RESULTS FROM FOCUS 

Stefano Bianco 
Laboratori Nazionali di Frascati dell'INFN 
on behalf 1 of the FOCUS Collaboration 



Abstract 

New results from the photoproduction experiment FOCUS are reported: Dalitz plot analysis, semileptonic 
form factor ratios and excited meson spectroscopy. 
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I report 1 on three new results from the photoproduction experiment FOCUS: the first Dalitz plot 
analysis of charm meson decays using the K-matrix approachf [Q, new measurements of the Df — > 
0(1020) n + v form factor ratios[|2], and new measurements on L=l excited meson spectroscopy[ 3 1, i.e., 
precise measurements of the masses and widths of the and D^ mesons, and evidence for broad states 
decaying to D + ir~ , D°tt + (the first such evidence in D°ir + ). The data for this paper were collected in the 
Wideband photoproduction experiment FOCUS during the Fermilab 1996-1997 fixed-target run. 

1 Dalitz plot analysis of Df and D + decay to tt + -k ~tt + using the K-matrix formalism 

Charm-meson decay dynamics has been extensively studied in the last decade. The analysis of the three- 
body final state by fitting Dalitz plots has proved to be a powerful tool for investigating effects of resonant 
substructure, interference patterns, and final state interactions in the charm sector [ 4 5 1. The isobar formal- 
ism, which has traditionally been applied to charm amplitude analyses, represents the decay amplitude as a 
sum of relativistic Breit-Wigner propagators multiplied by form factors plus a term describing the angular 
distribution of the two body decay of each intermediate state of a given spin. Many amplitude analyses 
require detailed knowledge of the light-meson sector. In the case of a narrow, isolated resonance, there is a 
close connection between the position of the pole on the unphysical sheet and the peak we observe in exper- 
iments at real values of the energy. However, when a resonance is broad and overlaps with other resonances, 
this connection is lost. The Breit-Wigner parameters measured on the real axis (mass and width) can be 
connected to the pole-positions in the complex energy plane only through models of analytic continuation. 

A formalism for studying overlapping and many channel resonances has been proposed long ago 
and is based on the K-matrix | (i 7 1 parametrization. This formalism, originating in the context of two-body 
scattering, can be generalized to cover the case of production of resonances in more complex reactions [ 8 1, 
with the assumption that the two-body system in the final state is an isolated one and that the two particles 
do not simultaneously interact with the rest of the final state in the production process [ 7|. The K-matrix 
approach allows us to include the positions of the poles in the complex plane directly in our analysis, thus 
directly incorporating the results from spectroscopy experiments. 

Full details on event selection and analysis cuts are reported in [[T). The Dalitz plot analyses are 
performed on events within 2 a the nominal Df or D + mass (Fig. The decay amplitude of the D 
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Figure 1: Df (left) and D + (right) Dalitz plots. 

meson into the three-pion final state is written as A(D) = a e lS " + Fi + J2i a,ie lSi B(abc\ri) where the 
first term represents the direct non-resonant three-body amplitude contribution, F\ is the contribution of 
S-wave states and the sum is over the contributions from the intermediate two-body non-scalar resonances. 
B(abc\ri) are Breit-Wigner terms. The amplitude for the particular channel (00)^ + tt can be written in the 
context of the K-matrix formalism as F\ = (I — iKp)jrPj where / is the identity matrix, K is the K-matrix 
describing the isoscalar S'-wave scattering process, p is the phase-space matrix for the five channels, and 
P is the "initial" production vector into the five channels. In this picture, the production process can be 
viewed (Figure[2) as consisting of an initial preparation of several states, which are then propagated by the 
(I — iKp)^ 1 term into the final one. Only the F\ amplitude is present in the isosinglet S'-wave term since 
we are describing the dipion channel. 

We use the K-matrix parametrization of (00) ++ -wave scattering following obtained through a global 
fit of the available scattering data from itit threshold up to 1900 MeV, see [ 9 1. The results are presented in 
Tabled 

In conclusion, the K-matrix formalism has been applied for the first time to the charm sector in our 
Dalitz plot analyses of the Df and D + — > ir + ir~ir + final states. Furthermore, the same model is able to 
reproduce features of the D + — » tt + tt^tt + Dalitz plot that otherwise would require an ad hoc a resonance. 
In addition, the non-resonant component of each decay seems to be described by known two-body S-wave 
dynamics without the need to include constant amplitude contributions. 

The K-matrix treatment of the S-wave component of the decay amplitude allows for a direct inter- 
pretation of the decay mechanism in terms of the five virtual channels considered: ttit, KK, rjrj, rjrj' and 
4 7r. The resulting picture, for both Df and D + decay, is that the S-wave decay is dominated by an initial 
production of r\r\, rjrj' and KK states. Dipion production is always much smaller. This suggests that in 
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Table 1: Results on and D + — > n + ir~Tt + fit fractions and phases. Beside the first reported error, which 
is statistical, two systematic errors are quoted. The first one is from the measurement systematics and the 
second one is due to the particular solution chosen for the K-matrix poles and backgrounds. 



D+ 


decay channel 


fit fraction (%) 


phase (deg) 


(S-wave) 


87.04 ±5.60 ±4.17 ±1.34 


(fixed) 


/ 2 (1270)tt+ 


9.74 ±4.49 ±2.63 ±1.32 


168.0 ±18. 7 ±2. 5 ±21. 7 


p°(1450)7r+ 


6.56 ±3.43 ±3.31 ±2.90 


234.9 ± 19.5 ± 13.3 ±24.9 


D+ 


decay channel 


fit fraction (%) 


phase (deg) 


(S-wave) 7r+ 


56.00 ±3.24 ±2.08 ±0.50 


(fixed) 


/2(1270)vr+ 


11.74 ±1.90 ±0.23 ±0.18 


-47.5 ±18.7 ±11.7 ±5.3 


P (770)tt+ 


30.82 ±3.14 ±2.29 ±0.17 


-139.4 ±16.5 ±9.9 ±5.0 




Figure 2: K-matrix picture of a D meson decay to three pions, with a dipion in a isosinglet S-wave. 



both cases the S-wave decay amplitude primarily arises from a ss contribution such as that produced by 
the Cabibbo favoured weak diagram for the and one of the two possible singly Cabibbo suppressed 
diagrams for the D + . For the D + , the ss contribution competes with a dd contribution. That the /o(980) 
appears as a peak in the tttt mass distribution in D + decay, as it does in Df decay, shows that for the 
S-wave component the ss contribution dominates. Comparing the relative S-wave fit fractions that we 
observe for and D + reinforces this picture. The S-wave decay fraction for the D+ (87 %) is larger 
than that for the D + (56 %). Rather than coupling to an S-wave dipion, the dd piece prefers to couple 
to a vector state like /3°(770) which accounts for ~ 30% of the D + decay. This interpretation also bears 
on the role of the annihilation diagram in the D+ — ► 7r + 7r~7r + decay. Our data suggest that the S-wave 
annihilation contribution is negligible over much of the dipion mass spectrum. It might be interesting to 
search for annihilation contributions in higher spin channels, such as p°(1450)7r and /2(1270)7r. 

2 New measurements of the D+ — > ^(1020) fi + is form factor ratios 

The Df — » 0(1020) n + v decay amplitude is described by four form factors with an assumed (pole form) 
q 2 dependence. The Df — » 0(1020) n + v amplitude is then described by ratios of form factors taken at q 2 
= 0. The traditional set is: r^, r%, and r v . According to flavor SU(3) symmetry, one expects that the form 
factor ratios describing Df — » 0(1020) should be similar to those describing D + — > K* (892)° n + u 
since the only difference is an s spectator quark instead of a d spectator quark. The existing lattice gauge 
calculations [ 10 1 predict that the form factor ratios describing £>+ — > 0(1020) t + V( should lie within 10% 
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Figure 3: Form-factor ratios comparison in previous data, and the new FOCUS result. World averages (not 
including this result) are also shown (shaded bands). 



of those describing D + — > K*(892)°£ + ve. Although the measured r v form factors are quite consistent 
between D+ -> 0(1020) and D+ -> K*(892)°£ + vj>, there is presently a 3.3 a discrepancy between 
the T2 values measured for these two processes with the previously measured Df — ► 0(1020) l + vg value 
being a factor of about 1.8 times larger than the r 2 value measured for D + — > A"*(892)°^ + t'f ( Figure[3}. 
For a review and references see [ 4|, while full details on event selection are found in [ 2\. The m K + K - 
distribution for the Df — > K + K~ p + v candidates is shown in Figure^] 

The r v and T2 form factors were fit to the probability density function described by four kine- 
matic variables (q 2 , cos#v, cosOe, and x) f° r decays in the mass range 1.010 < m K + K - < 1.030. 
We find r„=1.549 ± 0.250 ± 0.145, r 2 =0.713 ± 0.202 ± 0.266. Our measured r v and r 2 values for 



200 
^ 150 
>100 
50 


1.00 1.05 1.10 

M(K + K") (GeV/c 2 ) 

Figure 4: The data is the solid histogram and cc background Monte Carlo is the dashed histogram. The cc 
background Monte Carlo is normalized to the same number of events in the sideband region 1.04 GeV/c 2 
< mK+K- < 1-14 GeV/c 2 . 
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D+ — > 0(1020) [i^v are very consistent with our measured r v and values for D + — ► X*(892) /i + i/ [ 
1111 . The measurements reported here call into question the apparent inconsistency between r-2 values the 
D + K*(%92)°l + vi and £)+ — > 0(1020) form factors present in previously published data and 
are consistent with the theoretical expectation that the form factors for the two processes should be very 
similar. 

3 L=l excited charm meson spectroscopy 

High-statistics datasets from fixed-target and e + e~ colliders have recently provided the physics community 
with a wealth of data on excited charm meson spectroscopy. I report in this paper on new results [ 3 1 on L=l 
cu, cd, cs states, pointing the reader to detailed reviews for an account of the experimental scenario ri!2ll4l. 
In the limit of infinitely heavy quark mass, the heavy-light meson behaves analogously to the hydrogen 
atom, i.e., the heavier quark does not contribute to the orbital degrees of freedom (which are completely 
defined by the light quark). The angular momentum of the heavy quark is described by its spin Sq, and that 
of the light degrees of freedom are described by j q = s q + L, where s q is the light quark spin and L is the 
orbital angular momentum of the light quark. The quantum numbers Sq and j q are individually conserved. 
The quantum numbers of the excited L = 1 states are formed by combining Sq and j q . For L = 1 we have 
j q = 1/2 and j q = 3/2. When combined with Sq they provide two j q = 1/2 (J=0,1 where J is the total 
angular momentum of the excited charm meson) states, and two j q = 3/2 (J=l, 2) states. In this paper these 
four states will be denoted by D% , D i ( j q = 1 / 2 ) ,D i ( j q = 3/ 2 ) and D\ . 

Analysis procedures are explained in detail in [|3]. The L = 1 charm mesons were reconstructed 
via D + ir~ and D°ir + combinations. The D° decays were reconstructed in the channels D° — > K~ir + and 
D° — > K~Tr + ir + -K~ . The D + decays were reconstructed in the channel D + — > K~n + ir + . Our starting 
samples for these decay modes are 210,000, 125,000 and 200,000 events, respectively. Figure|5J;) shows 
the distribution of the invariant mass difference 

AM = M{(K--k + tt + )tt-) - M(K-tt + tt+) + M PDG (D+) 

where Mpdg(^ )+ ) is the world average D + mass [ 13|. Figure |5j;) shows a pronounced, narrow peak 
near a mass M ~ 2460 MeV/c 2 , which is consistent with the _D|° mass. The additional enhancement at 
M « 2300MeV/c 2 is consistent with feed-downs from the states and D^ decaying to D* + tt~ when 
the D* + subsequently decays to a D + and undetected neutrals. 
The mass difference 

AM+ = M((A"~7r + , K~tt + tt~tt + )tt + ) — M(K~ n + , K~ n + ir~ n + ) + Mpdg(-D°) 

spectrum (Figure|3Jl) shows similar structures to the AAfo spectrum. The prominent peak is consistent with 
a £>2 + of mass M ~ 2460 MeV/c 2 . The additional enhancement at M « 2300 MeV/c 2 is again consistent 
with feed-downs. 
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Table 2: Measured masses and widths for narrow and broad structures in D + ir and D°ir + invariant mass 
spectra. The first error listed is statistical and the second is systematic. Units for the masses and widths are 
MeV/c 2 . 







D*+ 


D° 




Yield 


5776 ± 869 ± 696 


3474 ± 670 ± 656 


9810 ± 2657 


18754 ± 2189 


Mass 


2464.5 ± 1.1 ± 1.9 


2467.6 ± 1.5 ± 0.76 


2407 ± 21 ± 35 


2403 ± 14 ± 35 


PDG03 


2458.9 ±2.0 


2459 ± 4 






Width 


38.7 ±5.3 ±2.9 


34.1 ±6.5 ±4.2 


240 ± 55 ± 59 


283 ± 24 ± 34 


PDG03 


23 ±5 


25± 8 7 







We fit the invariant mass difference histograms with terms for the D^ , peaks, D\ and D'^ 
feed-downs, combinatoric background and the possibility of a broad resonance. The broad resonance is 
necessary to obtain a fit to the data of acceptable quality (Fig[5]c-d). Our final results are shown in Table|2] 
Our mass measurement of the broad state is higher than a recent measurement by BELLE [ 14 1. Our result 
on the broad state have stimulated a series of theory studies, which try to reconcile the experimental picture 
of excited non-strange, and strange charmed mesons. 
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